Eddington's theory of gravity and its progeny. 
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We resurrect Eddington's proposal for the gravitational action in the presence of a cosmological 
constant and extend it to include matter fields. We show that the Newton-Poisson equation is 
modified in the presence of sources and that charged black holes show great similarities with those 
arising in Born-Infeld electrodynamics coupled to gravity. When we consider homogeneous and 
isotropic space-times we find that there is a minimum length (and maximum density) at early times, 
clearly pointing to an alternative theory of the Big Bang. We thus argue that the modern formulation 
of Eddington's theory, Born-Infeld gravity, presents us with a novel, non-singular description of the 
Universe. 



I. INTRODUCTION 

The Einstein-Hilbert action has been the mainstay of 
gravitational theory for almost a century. It can be ex- 
pressed as Seh = I / d A XyJ\g\(g al3 R a p — 2 A) where g a p 
is the metric of space-time (and \g\ is its determinant), 
R aj 3 is the Ricci tensor of that metric, A is the cosmolog- 
ical constant (in this letter we will work in Planck units 
8ttG = 1). An intriguing, alternative proposal for the 
gravitational action was proposed by Eddington in 1924 
Pi. He suggested that, at least in free, de-Sitter space, 
the fundamental field should be the connection, and 
the relevant action should be 



S 



Edd 



= 2k, 



(1) 



where K is a constant with inverse dimensions to that of 
A and \R\ is the determinant of R^. Note that R a p is 
constructed solely from the connection. (Here and in the 
rest of this Letter, R^ V {T) represents the symmetric part 
of the Ricci tensor build with the connection.) 

Varying Sem-i integrating by parts and eliminating a 
vanishing trace we obtain 



V q (2kVI- r I^) = 



to General Relativity. In fact they can be viewed as be- 
ing dual to each other- while Seh is proportional to A, 
^Edd is inversely proportional to A and one can imagine 
that they should be useful in different regimes. 

Eddington's theory of gravity is incomplete in that it 
doesn't include matter. There have been subsequent at- 
tempts at coupling matter to Y in the original Eddington 
spirit 0]. The idea is as follows. Start with a Palatini 
gravitational action coupled to matter I[g,Y,^>). Here 
$ denotes all matter fields. The metric enters with no 
derivatives and it is, in principle Q, possible to express 
g^u in terms of Y and $ using its own equation of mo- 
tion. Since this step is algebraic, one can replace this 
formula back into the action and obtain an "affine varia- 
tional principle" i 7 [r,\t'] depending only on the connec- 
tion and matter fields. The action /'[r, 1 !/] derived by 
this procedure may be quite complicated |2| , but it must 
be stressed that its dynamics is fully equivalent to the 
original metric theory. As a matter of fact, this proce- 
dure is most simple without matter and precisely maps 
the Palatini version of Seh into SEdd- 

In this Letter, we shall reconsider the problem of cou- 
pling Eddington action to matter without insisting nei- 
ther on a purely affine action nor on a theory equiva- 
lent to Einstein gravity. We shall allow the metric to 
be present and switch to a Born-Infeld [H like structure, 



where V is the covariant derivative defined in terms of v / lsV + >^V|. Observe that for large kR^ one reob- 



rjg 7 and is the inverse of the Ricci tensor. This equa- 
tion can be partially solved if we define a new rank-2 ten- 
sor such that V Q (-\/|g|(7 A1!y ) = 0. The field equations 
for this theory then become 



tains the Eddington functional. A purely metric proposal 
can be found in [5[ . We shall focus on a Palatini formula- 
tion (with g^ w and Y^ a/3 independent) of the gravitational 
action first proposed in 



2KyJ\R\R^ 
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(2) S BI [g,Y,*} 



d 4 x 



\J\g^ v + kR^AY)\ - X^/g 



which can be rewritten as the Einstein field equations if 
we equate q a p with g a p and k with A . Hence, Edding- 
ton's action should be a viable, alternative starting point 



s M \g,r,v] 



(3) 
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where A is dimensionless. Note that A must be different 
from zero. For A = the metric variation yields (with 
no matter) y/\g + nR\ [g + kR] 1{iV = which clearly 
makes no sense. In @ the matter fields were introduced 
in a non-conventional way inside the square root. Here 



2 



we add matter in the usual way. 

This action has all the correct limits and no patholo- 
gies. For small values of kR, the action ((3]) reproduces 
the Einstein- Hilbert action with A = (A — 1)/k. On 
the other hand for large values of nR, the action ap- 
proximates Eddington's. Unlike the matter actions built 
in Q and reviewed above, ([3]) is not equivalent to the 
Einstein-Hilbert action and we shall observe interesting 
deviations especially for large curvatures when Edding- 
ton action dominates. 

In vacuum with Sm = 0, on the other hand, the action 
([3]) is equivalent to Seh- This can be proven in two steps. 
With Sm = the equation of motion for g^ v implies 
9iiv = xrr-fiW- We then replace the metric back into the 
action to obtain the Eddington action -^y J y/\R\ which 
we already remarked is equivalent to Seh- The case A = 
1 is not degenerate although this particular proof fails. 
For A = 1 the equivalence can easily be proven at the 
level of the equations of motion. This is to be contrasted 
with f(R) theories Q, either metric or Palatini, where 
deviations are observed even in vacuum. We shall now 
consider ([3]) with matter in different situations to explore 
it properties and we shall restrict ourselves to the case 
where matter only couples to the metric (couplings to T 
may arise due to quantum gravitational corrections). See 
7] for other discussions. 

The equations of motion for this theory are the follow- 
ing. Varying with respect to g^ v one obtains, 



V\9 + kR \ 



[(g + kR) 



(4) 



Here T^ v is the standard energy momentum tensor with 
indices raised with the metric g^. 

The variation with respect to T can be simplified by 
introducing an auxiliary metric q^ v compatible with V. 
The equation of motion becomes 



and (|4|) one finds the equation 

>/g«T = X V99^ - Ky/gT'"'. 



(5) 

9ap,<r)- Combining (O 
(6) 



Again, here we have that q^ u is the inverse of 
q^ v . Equations (0 and |6]) form a complete set and pro- 
vide the simplest set of equations to study this theory. 
The conservation equations for the matter fields are the 
same as in GR, = 0, where the covariant derivative 
here refers to the metric g^ v . Although this conserva- 
tion equation is not at all obvious from the equations 
of motion it is nevertheless true because matter is 
covariantly coupled to g^. 

We can now rewrite equation (|4| as 
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kRI = 



«A)g- 1 - Kg - 1 T g - 1 ][ g - 



cRl 



where, for the moment, we are writing the equation in 
the matrix format, I is 4 x 4 identity matrix. If we take 



the determinant of both sides and replace it in the field 
equation, we find 

g + kR = I^Kl + KAJg^-Kg^Tg- 1 !* 

xKl + KAJg^-Kg^Tg- 1 ]" 1 

i.e. we can solve for R^ v in terms of g^y and T^ v . Note 
that this equation is generic. Recall however that is 
a function of T° (and therefore q^ v ). 

Expanding the field equations to 2 nd order in n to find 
the 1 st order corrections to Einstein's equations 



R 



S. 



1 



:Sgnv 



where S^ = T^T av — \TT^ V (note that A does not con- 
tribute to S). Combining this equation with the second 
field equation, ([5]), we have the lowest order correction 
to Einstein gravity. 

With the lowest order correction in hand, we can study 
the non-relativistic limit. As expected, the Poisson equa- 
tion is modified in this theory. Consider a time indepen- 
dent metric ds 2 = -(1 + 2$)di 2 + (1 - 2^)dx ■ dx where 
$ and $ only depend on x, and an energy momentum 
tensor T^ u = pu^u v . We linearize equations ((H) and ([5]) 
keeping terms linear in $, "J" and p. The full set of lin- 
earized equations are solved by $ = \P plus the modified 
Poisson equation, 



1 



1 



(7) 



V 2 $ = — p kV 2 a. 

2 F 4 H 

We come up against a key characteristic of this theory- it 
reproduces Einstein gravity precisely within the vacuum 
but deviates from it in the presence of sources. We expect 
it therefore to play a role in regions of high density, such 
as within a black hole or in the very early universe. 

An obvious next step is to explore the inner structure 
of black holes. Let us consider spherically symmetric 
configurations of the form 

dr 2 

ds 2 = -i){rff{r)dt 2 + — + rW. (8) 
f(r) 

As we have just remarked, in vacuum, the action (J3J 
is fully equivalent to Seh with A = (A — 1)/k. The 
Schwarzschild-dS metric with tp = 1 and / = 1 — 2m/ r — 
Ar 2 /3 is thus a solution to ([3]) with no sources. New 
interesting effects show up when matter is present. 

The simplest form of "matter" is an electromagnetic 
field (matter in the sense that ^ 0). Let us add 
— \\fg~F ilu F liV to the action (|3j). The exact solution can 
be found and involves elliptic integrals. For simplicity we 
set A = 1 yielding an asymptotically flat geometry. The 
solution is: 



m 

ip(r) 
E(r) 



ry /, r i + Kq 2 
— nq z 



^2r 4 + 2nq 2 
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\/r 4 + Kq 2 



dr 



(r 2 — q 2 )(r 4 — nq ) 
r 4 y/r 4 + Kq 2 



2M 
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where E(r) is the electric field. M is the mass and q the 
electric charge. Observe the similitude with pure Born- 
Infeld electrodynamics. For k > the electric field is 
everywhere regular. The metric functions and still sin- 
gular at r = and at r 2 = y/nq. This singularity is 
however inside the horizon. 

Strictly speaking, this solution represents the exterior 
solution to a charged object. In this sense, the r = 
singularity does not yet imply a singular solution. In 
order to explore the singularity structure of the action ^ 
one needs to add normal matter and study a collapsing 
object. A simple model is an interior cosmology glued 
together with an exterior Schwarzschild metric. We shall 
discuss the details of this problem elsewhere but we can 
anticipate interesting conclusions: as we now show the 
cosmology associated to ([3]) predicts a maximum density 
and a singularity free Universe. 

It has become clear that Eddington gravity can play a 
role in regions of high curvature or density, the conditions 
one might expect in the early Universe. We therefore 
focus now on cosmology. We shall assume a homogeneous 
and isotropic metric, 



ds 2 



-dt 2 + a(t) 2 dx ■ dx 



coupled to an ideal fluid T>* v = (p + p)u^u u + pg 1 *" '. As 
we remarked before, the fluid satisfies its standard con- 
servation equation p — —3H(P + p) 
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FIG. 1: The Hubble rate H 2 as a function of energy density, 
p, for a radiation filled Universe with no cosmological con- 
stant and |re| = 1 We normalize the density by pb (for which 
H 2 (pb) = 0); for re > we have ps = 3/re while for re < we 
have pB = — 1/k. 

We can assume that goo = —U and qij — a 2 VSij and 
use equations (|U and flSJ) to find 

D , „ D 



U = 



1 



and V = 



np T 



1 - kP t 



where D = + npr){l — kPt) 3 , Pt = p+A and Pt 
P — A. If we assume that P = wp, we can define 



f(pA) 



1 



3{k Pt + kPt)(1 



w 



npr — kPt) 



(l + K.p T )(l -kP t ) 



G(p,A) = - 

K 



-4 




FIG. 2: The scale factor, a, normalized by the minimum 
length, aB, as a function of time, t (in arbitrary units). In 
the top panel we show the existence of a bounce for re < 
while in the bottom panel, the scale factor approaches the 
minimum length as a — > — oo. 



to find the Friedman equations: 

8G_ 

9F 2 



(9) 



If we assume w — (i.e. a dust filled universe with 
a cosmological constant), we can expand equation ^ in 
terms of reA to find corrections to the late time Friedman 
equations 



3H Z 



(P + A) + 



T-<" + A > 



KA + 0[(reA) s 



As expected, we find that we recover conventional Fried- 
man cosmology at late times. 

Let us now focus on the evolution of the scale factor 
at early times. Assuming radiation domination, we have 
P t = p, Pt = P = \p and, definining p — up we find 
that equation © becomes 



m\p)= - 



i 



:: (i + p)(3-p) 2 

(3 + p 2 ) 2 

For small p we recover the conventional Friedman uni- 
verse, H 2 — p/3 but at high densities we come up against 
a novel effect: we find a stationary point, H 2 = at p = 3 
(for k > 0) and at p = 1 (for re < 0). The new stationary 
points correspond to a maximum density, pg; in Figure 
(JTJ we plot the H 2 as a function of p within the physically 
acceptable region. 

If k is of order unity, then we can interpret this effect 
as a cutoff in energy density at around the Planck scale. 
Given that p oc a~ 4 (recall that the conservation equation 
is the same as in GR, so this relation does not change in 
this theory) , this means that there is a minimum value for 
the scale factor at a B ~ 10 _32 (re) 1 / 4 a (where a is the 
scale factor today)- corresponding to a minimum length, 
aB = (/?o//3b)~ 4 , in cosmology. 
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The nature of the expansion rate at as depends on the 
sign of n as we show in Figure [5] where we plot the scale 
factor, a, for both situations. For k < 0, one can show 
that H 2 oc (a— as) which means that a— as oc \t— ts\ 2 for 
a fixed ig. This means that at is the Universe undergoes 
a regular bounce. The transition through the bounce 
occurs on a time scale of At ~ y/H. Again, if k is if order 
unity, the deviation from conventional cosmology occurs 
at the Planck time on a time scale of that order. 

A far more interesting, in our view, behavior can be 
found if K > where one can show that H 2 ~ (a — as) 2 
which means that ln(a/as — 1) oc f — ts- In this case 
there is no bounce; if we wind back the clock, the en- 
ergy density will reach a point (corresponding to about 
ps/2 as can be seen from Figure [T]) in which acceler- 
ated expansion kicks in. As above, this corresponds to 
what we would perceive as the Planck time, but now, 
the scale factor will take an infinitely long time to reach 
as- This leads to an intriguing alternative cosmology at 
early times in which, depending on how close its initial 
density is to pg, the Universe loiters for a long time (in 
terms of Planck units) until it emerges into a standard 
cosmological evolution. This evolution can be seen as 
an alternative incarnation of the Einstein static universe 
at early times. As in the conventional construction, the 
loitering phase is unstable to expansion. 

There are two effects that may be of importance in 
the early Universe. First of all, it should be clear that we 
have been looking at the classical behaviour of this theory 
of gravity. The onset of a bounce or minimum length 
may signal pathologies at the quantum level, such as the 



presence of ghosts and negative normed states. A more 
detailed analysis of the gravitational Born-Infeld action 
will allow us to check if these pathologies indeed exist. 
Secondly, we have discarded the effect of A (i.e. kA <C 1) 
in our analysis of the Early Universe. If, however, there 
is an early time contribution to A from, for example, 
an infiaton, the dynamics will be different; with kA ~ 1 
there will be another inflection point and H 2 (p) will never 
reach for a finite density. Such scenario merits further 
analysis. 

Eddington's theory of gravity, through its descendant, 
the gravitational Born-Infeld action, may lead to an en- 
tirely new view of the Universe. As pointed out above, 
the presence of a maximum density in cosmology may 
have a bearing not only on the early universe but in the 
dynamical formation of black holes. Thus, if a mimimum 
length also arises during gravitational collapse, the Uni- 
verse may be entirely singularity free [i~l| , solving one of 
the problems that has troubled relativists since Einstein 
first proposed his theory of gravity. 
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